MeV-mass dark matter and primordial nucleosynthesis 
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The annihilation of new dark matter candidates with masses mx in the MeV range may account 
for the galactic positrons that are required to explain the 511 keV 7-ray flux from the galactic bulge. 
We study the impact of MeV-mass thermal relic particles on the primordial synthesis of ^H, *He, 
and ^Li. If the new particles are in thermal equilibrium with neutrinos during the nucleosynthesis 
epoch they increase the helium mass fraction for mx ^ 10 MeV and are thus disfavored. If they 
couple primarily to the electromagnetic plasma they can have the opposite effect of lowering both 
helium and deuterium. For mx = 4-10 MeV they can even improve the overall agreement between 
the predicted and observed and *He abundances. 

PACS numbers: 95.35.+d, 26.35.+C, 14.80.-j 



I. INTRODUCTION 

In a recent series of papers the intriguing possibil- 
ity was explored that the cosmic dark matter consists 
of new elementary particles with masses in the MeV 
range [H 13, 13, S 111 ■ While weakly interacting massive 
particles (WIMPs) as dark matter candidates, notably 
supersymmetric particles, are usually thought to have 
masses exceeding tens of GeV, it is not difficult to come 
up with viable MeV-mass candidates such as sterile neu- 
trinos or axinos |7l|. The remarkable aspect of the 
MeV candidates studied in Refs. 0, S S |1 is that 
they are taken to be thermal relics and thus require a 
primordial annihilation rate much larger than given by 
ordinary weak interactions. Contrary to naive intuition, 
such particles are not excluded by any obvious laboratory 
measurement or astrophysical argument [E S 13, S IS 13 ■ 
Quite on the contrary, it was argued that the annihilation 
of these dark matter particles in the galactic bulge can 
produce enough positrons to explain the 511 keV 7-ray 
signature that was recently confirmed by the INTEGRAL 
satellite and that seems difficult to explain with 

traditional astrophysical sources. Alternatively, decaying 
low-mass dark matter particles have also been proposed 
as a positron source [13, llll • 

In the early universe, weak interactions freeze out at 
a temperature of about 1 MeV, just before the epoch of 
big-bang nucleosynthesis (BBN). Therefore, MeV-mass 
particles with larger-than-weak interaction rates are ex- 
pected to be in thermal equilibrium at BBN and thus 
would add to the primordial energy density and expan- 
sion rate. One naively expects such particles to increase 
the primordial helium mass fraction Yp , exacerbating the 
tension between observations and the BBN prediction. 
This expectation bears out for new particles coupling 
primarily to neutrinos. However, the new particles pro- 
posed in Refs. [1 S H S would primarily interact 
with the electromagnetic plasma. It turns out that the 
BBN effects of these particles are non-trivial and cer- 
tainly disfavor masses below about 2 MeV, but masses 
in the approximate range 4-10 MeV actually have the 
opposite effect of lowering Yp without significantly af- 



fecting deuterium. Therefore, the concordance between 
the predicted and observed '*He and abundances is 
slightly improved. This little-known but intriguing effect 
was found a long time ago in a general study of the im- 
pact of new particles on BBN [I^- The main purpose 
of our work is to re-examine this effect in the context of 
the MeV-mass dark matter hypothesis and with the help 
of a modern BBN calculation and current observational 
data. 

In addition to our BBN study presented in Sec. II we 
briefly consider two other possible consequences of MeV- 
mass dark matter particles. In Sec. Ill we note that the 
same mechanism proposed to explain the 7-ray signa- 
ture from the galactic bulge should also produce a dif- 
fuse background of low-energy cosmic-ray positrons in 
the solar neighborhood. This argument strengthens the 
conclusion reached independently in Refs. ^ 2tc ^ ^ ^ 
that an s-wave annihilation cross section for the reaction 
XX e+e~ is strongly disfavored. Further, in Sec. IV 
we consider possible energy-transfer effects caused by 
such particles trapped in the Sun or other stars. These 
effects turn out to be very small because the trapping 
efficiency is inhibited by the smallness of mx- In Sec. V 
we summarize our conclusions. 



II. MEV-MASS PARTICLES AND BIG-BANG 
NUCLEOSYNTHESIS 

A. Impact of New Particles 

Thermal relic particles freeze out at a cosmic temper- 
ature Tp that for weakly interacting particles is about 
1 MeV. If the particle mass is somewhat larger than Tp 
the number density is suppressed by annihilations so that 
the relic density is maximal for mx of order Tp, i.e. for 
neutrinos the relic mass density would be maximal for 
MeV-range masses and overdose the universe by about 
five orders of magnitude ^^l- Therefore, MeV-mass ther- 
mal relics must interact far more strongly than neutrinos 
so that they annihilate more efficiently to reduce their 
density to a level compatible with the dark-matter abun- 
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dance. We briefly summarize the relationship between 
the annihilation cross section and the relic density in Ap- 
pendixEl The result reveals that MeV-mass dark matter 
particles must have been in thermal equilibrium through- 
out most of the primordial nucleosynthesis epoch so that 
the exact final dark-matter abundance is not relevant for 
our study. We will use the equilibrium assumption in our 
numerical implementation that was performed by a mod- 
ified version of the new BBNCODE recently developed in 
Naples and documented in Ref. . 

The impact of our new particles on BBN differs qual- 
itatively and quantitatively depending on the dominant 
annihilation and scattering channels. Before turning to 
specific cases we discuss how the BBN equations and in- 
puts are affected. To this end we first note that the 
new particles enter BBN through the second Friedmann 
equation H'^ — (Stt/S) Gn Ptot by their contribution to 
the total mass-energy density ptot. We consider the X- 
particles to be in perfect thermal equilibrium and assume 
that the number densities of particles and antiparticles 
are the same for those cases where X ^ X. The number 
density, pressure, and energy density contributed by the 
new particles is thus taken to be 
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where x = mx/Tx and the sign + (— ) refers to fermion 
(boson) statistics. Tx is identical with the ambient neu- 
trino temperature T^, for X-particles that dominantly 
couple to neutrinos, while Tx = T-y for the electromag- 
netically coupled case. Since y > x > and 
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we may expand the thermal integrals as 
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where Ki{x) is the special Bessel function of order I. In 
the numerical code the series are truncated at n = 4. 
The error of the physical quantities is always < 0.18%. 

A more subtle effect is caused by the conservation of 
entropy that implies a modified relation between the neu- 
trino temperature and that of the electromagnetic 
plasma T = T^. The cosmic entropy density is 
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The second equality applies for relativistic species, i.e. in 
the absence of warm species where m ~ T. Note that 
gl^ is always considered as a function of the temperature 
T of the electromagneticplasma, not of Ti (a different 
notation is used in Ref. [131 )■ After neutrino decoupling 
at Td ~ 2.3 MeV the entropy in a comoving volume 
is separately conserved for the "neutrino plasma" and the 
electromagnetic one so that 
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For new particles being only electromagnetically coupled 
this equation simplifies to 
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thus giving a higher temperature ratio relative to the 
standard case. This simulates the effect of N^s < 3 ther- 
mally excited neutrino species at BBN and thus a reduced 
primordial helium abundance. The opposite is true for 
j'-coupled new particles. 

The modified T^{T) relation also affects the n ^ p 
weak rates that depend on both T and T^, through the 
phase-space dependence of the initial and final states of 
the processes. We implemented the modification of these 
rates in a perturbative way by introducing the small pa- 
rameter 
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where T^{T) is the standard dependence on the electro- 
magnetic temperature T. Typically 6 assumes values of 
order 0.01 and is always smaller than about 0.1. The neu- 
trino temperature enters the weak rates through Fermi 
factors of the kind 



1 

exp(azi,) 



(9) 



3 



where Zy = me/Ty Therefore, the additional terms for 
the rates can be obtained by integrating the factors 
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exp[azO(l + (5)] 1 + exp(az5J) 



(10) 

where in our numerical treatment the series was trun- 
cated to the third term. The corrections and the standard 
rates in the Born approximation were used to compute 
numerically the relative changes to the rates. Therefore, 
in the limit (5 — > one recovers the n <-> p rates including 
finite mass, QED radiative and thermodynamic correc- 
tions, while disregarding modifications of these sublead- 
ing effects in the 5^0 corrections. We finish with six 
functions ti{T) 
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fitting the change in the weak rates with an accuracy 
better than 5%. 

In principle, one has a single covariant energy con- 
servation equation for all components of the primordial 
plasma. For the sake of simplicity, however, in the previ- 
ous considerations the "two fluids entropy conservation" 
was used to obtain the T^(T) relation. We can now de- 
rive the evolution of the thermodynamical quantities by 
applying the covariant energy conservation law to one of 
the two plasmas, e.g. the electromagnetic one, so that 
the first Friedmann equation is 
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where H depends on ptot through the second Friedmann 
equation. If the additional species X couples to the elec- 
tromagnetic fluid, the T-i-relation is further affected by a 
modified (pcm + Pcm) factor, at least until the scattering 
freeze-out is reached. This has been roughly estimated 
to happen at T ~ 35 keV (Appendix IB|| . In the numer- 
ical code, the X-particles were considered to decouple 
instantaneously from the electromagnetic component of 
the plasma for T < 35 keV. Relaxing this assumption 
our results remain essentially unchanged because for such 
low temperatures nucleosynthesis has almost completely 
stopped and the X-particles have a negligible impact, at 
least for the interesting mass range. 

Another input parameter for the BBN calculation is 
the radiation density contributed by ordinary neutrinos 
which we fix to N^s = 3. 

Finally, we need the cosmic baryon density. The best- 
fit value from the temperature fluctuations of the cosmic 
microwave radiation as measured by the WMAP satellite 
is ^Bh'^ = 0.024 ±0.001 Including large-scale struc- 
ture data from the 2dF galaxy redshift survey shifts this 
result to rie/i^ = 0.023 ± 0.001, and including Lyman-a 



data further shifts it to 0.0226 ±0.0008 In our study 
we always use a fixed value of 



QBh-" = 0.023. 



(13) 



We have checked that for 2a variations of JIb/i^ our con- 
clusions remain essentially unchanged. Note that in our 
BBN code the final value of 77 = n-Q/n^ or of Q,-Qh^ is 
used to work out the one that enters the initial condition 
of the problem. We have 
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so that entropy conservation implies the well-known fac- 
tor 11/4 in the standard case. In general it is a factor 
depending on the X-particle properties and was numeri- 
cally evaluated. 



B. Neutrino-Coupled Particles 

As a first generic type of X-particles we consider 
particles that annihilate predominantly into neutrinos 
XX ^ vv. We explicitly study three cases, that of Ma- 
jorana fermions with a total of gx = 2 inner degrees 
of freedom (case F2), self-conjugate scalar bosons with 
gx = 1 (case Bl), and scalar bosons with a particle and 
anti-particle degree of freedom {gx = 2, case B2). With 
the ingredients discussed in the previous section we calcu- 
lated the abundances for the light elements ^H, ^He, and 
^Li shown in Fig. ^ as a function of the new particle mass 
mx- For mx ^ 20 MeV we recover the standard BBN 
predictions. For very small masses mx — > these par- 
ticles freeze out relativistically and their effect on BBN 
is exactly that of AiVcff = 4/7, 1, or 8/7 additional rel- 
ativistic neutrinos for the three cases Bl, F2 and B2, as 
summarized in Table 

The BBN effect of the new particles is dominated by 
their contribution to the primordial energy density and 
thus similar to an additional neutrino species. However, 
it is worthwhile to note the extra effect for intermediate 
mx relative to the asymptotic case mx — > 0. It is caused 
by the modified T^/T ratio previously described. 

Our results essentially agree, both qualitatively and 
quantitatively, with those of Ref. jl3l |. except for lithium. 
The difference is explained by our value of Q.^h'^ where 
the dependence of ^Li on A^oft is opposite from the situ- 
ation in Ref. I'^ o^'' "^^se this nuclide is essentially 
produced through the channel ^He(^He, 7)^Be(e~, i^e)^Li 



TABLE I: Cases for new particles 



Case 
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AiVeff 


Vi/Vf 


Bl 


Self-conjugate scalar boson 


1 


4/7 


3.25 


B2 


Scalar boson X ^ X 


2 


8/7 


3.75 


F2 


Majorana fermion 


2 


1 


3.65 
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0.26 



0.25 




N,ff=4 



Neff=3.5 
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TABLE 11: Measured primordial light-element abundances 
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FIG. 1: Calculated light-element abundances for *He (top), 
(middle), and ^Li (bottom) for neutrino-coupled new par- 
ticles. The indicated cases Bl, B2, and F2 are described in 
Table H] The horizontal lines indicate the effect of A'cff = 3, 
3.5, and 4 relativistic neutrinos. 



while at the lower value of flBh^ used in Ref. T3| the di- 
rect channel ^He(^H, 7)^Li dominates. 

Our theoretical predictions can be compared with the 
measured primordial abundances summarized in TablelTll 
For helium, the standard BBN prediction significantly ex- 
ceeds the most recent measured value To^ , and this dis- 
crepancy is even worse for other helium determinations 
that are lower (for a review see, e.g. Ref. HI). There- 
fore, while different observations of the primordial "^He 
abundance disagree on its exact value, a tension with the 
BBN prediction always exists. Our new particles exac- 
erbate this discrepancy for mx ^ 10 MeV and are thus 
disfavored or even excluded. 

The deuterium abundance extracted from the QSO 
systems agrees perfectly with the BBN prediction, al- 
though it could be affected by possibly underestimated 
systematic errors. In any event, the new particles change 
the prediction only within the la observational range so 



Element 



Abundance 



Helium" "^He Yp = 0.2421 ± 0.0021 



Reference 
[19] 

Deuterium D/H = 2.78+o;38 x 10"^ [2^] 
Lithium ^Li ''Li/H = (1.70 ± 0.17) x 10"^° Refs. in [15] 

Tor other determinations and a short review, see also ,18ij . 



that deuterium adds little new information on the viabil- 
ity of the new particles. 

A clear interpretation of the Spite plateau in the 
lithium data from metal-poor halo stars is still lacking 
so that it is not clear how to compare the lithium obser- 
vations (average value given in Table ^ with the the- 
oretical prediction that is roughly a factor 2-3 larger. 
Both standard and non-standard physics explanations of 
this discrepancy have been invoked, e.g. Ref. and 
references. Therefore, a meaningful comparison of our 
non-standard BBN prediction for lithium with observa- 
tions is difficult. We show these results primarily for the 
purpose of illustration. 



C. Electromagnetic Couplings 

We next turn to new particles that do not interact 
with neutrinos but remain in perfect thermal equilibrium 
with the electromagnetic plasma throughout the BBN 
epoch by virtue of processes such as XX ^ e~^e~ and 
X + e ^ X + e. The BBN impact of these particles is 
more subtle and can be opposite to the neutrino-coupled 
case. Moreover, this is the case relevant as a source for 
galactic positrons For particles of this sort, 

there is a tension between the required primordial anni- 
hilation cross section to achieve the correct dark-matter 
density and the one in the galaxy to avoid overproducing 
positrons. The solution prefered in Refs. .4, 5] is 

that of a predominant p-wave annihilation channel that 
will suppress the annihilation rate in the galaxy relative 
to that in the early universe. In particular, a specific 
model for a new boson was constructed where particles 
and anti-particles are not identical, i.e. our case B2. Of 
course, the detailed form of the cross section is not rele- 
vant for our work because we assume that the new par- 
ticles are in perfect equilibrium with the electromagnetic 
plasma. 

The calculated light-element abundances as functions 
of mx are shown in Fig. |21 For mx ^ 2 MeV the abun- 
dances are all shifted away from the observed values. 
However, for mx ^ 2 MeV the "entropy effect" works 
in the direction of lowering Yp, by up to AYp « —0.002 
for the B2 case, without significantly affecting the and 
^Li predictions. Therefore, the concordance with the ''He 
observations is improved. This point is illustrated more 
directly by Fig. |31 where we compare the new ^He predic- 
tions for the cases Bl and B2 with the standard N^s = 3 
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FIG. 3: ''He abundance as in Fig. El here on a linear scale 
for mx- The horizontal dotted line indicates the standard 
prediction for A'^off ~ 3. The gra^ band is the la observational 
range for Yp according to Ref. 
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We have always neglected the role of dark matter resid- 
ual annihilations during the freeze-out epoch. A detailed 
treatment of such effects is beyond the scope of our pa- 
per. However, an approximate study of this phenomenon 
(Appendix O indicates that this late-time entropy gen- 
eration effect is sub-dominant. Its effect goes in the di- 
rection of a further marginal reduction of Yp and a small 
increase in the deuterium yield. 

We have also neglected the possible photo-dissociation 
of and ^Li induced by late dark matter annihilations. 
As shown in Appendix IdI this effect is marginal unless, 
perhaps, if the dark matter particles couple directly with 
photons. 



FIG. 2: Calculated light-element abundances for ''He (top), 
^H (middle), and ^Li (bottom) for electromagnetically cou- 
pled new particles. The indicated cases Bl, B2, and F2 are 
described in Table The horizontal dotted lines indicate the 
standard predictions for the indicated values of the effective 
number of relativistic neutrino species, A'oft . 



case and with the recent observational determination [l9j 
for which the Icr error band is shown. Yp would benefit 
from this effect up to mx ^ 15 MeV, even if the val- 
ues mx ^ 10 MeV are prefered. Our results again agree 
qualitatively with those of Ref. ^l3| . 

For low masses, the value of rji needed to match the 
WMAP finding for r/f is significantly increased relative 
to the standard factor 2.75, as shown in the last column 
of Table ^ This explains physically the huge decrease 
(increase) in the ("^Li) yield that does not strongly 
depend on Nes- Quite on the contrary, as Yp depends 
only logarithmically on 77, its change is essentially domi- 
nated by the addition of extra degrees of freedom to the 
electromagnetic plasma. The V-particles are now hotter 
than in the neutrino case and thus provoke a bigger effect. 
For mx 20 MeV one recovers the standard predictions. 



D. Particles Coupled both to the Electromagnetic 
Plasma and to Neutrinos 

It may also be that the new particles interact strongly 
enough with both charged leptons and neutrinos to keep 
the two fluids in equilibrium beyond the usual decou- 
pling epoch, a situation that was not previously treated. 
In this case = is maintained longer and perhaps 
throughout the nucleosynthesis epoch, depending on the 
new particle properties. This effect would be present even 
for very high mx if one introduces an additional direct 
coupling of the neutrinos to the charged leptons through 
an extra gauge boson as in the prefered [/-boson model 
discussed in Refs. 0,0>ISSI3- Of course, in this case 
laboratory data, e.g. from v-e scattering j22, l23|, pro- 
vide strong limits so that this situation may be rather 
unphysical. 

In Fig. ^ we show the light-element abundances as a 
function of the assumed neutrino decoupling temperature 
Td. If it is lowered to values Tjy ^ me, the resulting mod- 
ification of the light-element abundances is quite extreme 
and strictly excluded. For and ^Li the "eta-effect" 
dominates, while for ^He the change in T^/T and that of 
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FIG. 4: Calculated light-element abundances for *He (top), 
(middle), and ^Li (bottom) if the neutrino and electro- 
magnetic plasmas are thermally coupled by a new interaction 
channel down to a temperature Tb. The horizontal lines in- 
dicate the standard predictions for the indicated values of the 
effective number of relativistic neutrino species, A'oft . 



III. LOW-ENERGY COSMIC RAY POSITRONS 

The proposed MeV-mass dark matter particles dis- 
cussed in Refs. 0, E S S are supposed to annihi- 
late in the galactic bulge and produce a flux of low- 
energy cosmic-ray positrons that can explain the ob- 
served 511 keV 7-ray signature. For such X-particles 
we also expect a flux of low-energy positrons at Earth 
from local dark-matter annihilation. Such a signature 
was already proposed for the detection of the traditional 
GeV-TeV mass range of dark matter particle candi- 
dates 0, For the case of annihilating MeV-mass 
dark matter particles we expect a much larger positron 
flux and therefore we study if additional constraints arise 
from the low-energy cosmic-ray positrons in the solar 
neighborhood. 

To this end we derive the expected positron flux from 
X-particle annihilation. For a stationary situation the 
continuity equation for the cosmic-ray positrons is 



n{E) d[b{E)n{E)] 



r{E) 



dE 



0, 



(15) 



where n{E) is the differential positron density, q{E) the 
injection spectrum, t{E) an effective containment time, 
and b{E) = —dE/dt the energy-loss function. This equa- 
tion is solved by 



i{E)^ 
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ds q{s) exp 



dy \ 
'{y)Kv)) 
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The expected positron flux at the top of the atmosphere, 
neglecting solar modulation effects, is je+ = n(£')/47r, 
where we have assumed relativistic velocities. 

In the relevant energy range and for an essentially neu- 
tral environment, ionization and bremsstrahlung are the 
dominant energy-loss mechanisms so that we have j26j 



h{E) « 36.18 meCTT^H 



H- 0.146 In 
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^ 709 me _ 



3.69 X 10"^^ MeV s^^nn [...]. 
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Pi, play a major role. 

If the neutrinos and electromagnetic fluids are only 
indirectly coupled by their interactions with the new 
X-particles, the limiting condition to prolong the ly- 
equilibrium with the photon bath is the number den- 
sity factor of the dark matter, nx, in addition to the 
interaction cross sections. Qualitatively, we expect that 



for mx 



< 



the combined effects conspire to hugely 



increase the ^He yield, while some compensation could 
appear in the and ''Li behavior. Obviously, for 
mx ^ 20 MeV one should recover the standard predic- 
tions, unless a direct v-e coupling exists. A detailed ki- 
netic study of this case would require a concrete particle 
physics model and is beyond the scope of our work. 



Here, ut is the Thomson cross section and tt-h the hydro- 
gen density that dominates the interstellar medium. For 
the relevant energies, b{E) is a slowly varying function. 

It is easy to estimate that low-energy positrons do not 
travel far before annihilating so that the containment 
time t{E) is identical with the annihilation time scale. 
The latter can be written in the form 123 



1.33 X 10^-^ s 



(18) 



where f{E) = 0.02-1 in the energy range of interest. 

Finally we need the positron injection spectrum from 
dark-matter annihilation. 



q{E)^nl{aav)S{l-E/mx).. 



(19) 
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where we have assumed that the X-particles are different 
from their antiparticles and that nx = nx- 

Equation H16|l cannot be expressed in closed form, but 
for our purposes an analytical approximation is accurate 
enough. The factor br does not depend on nn and is 
found to be 



br ^ 50-5000 MeV , 



(20) 



where the range reflects the monotonic energy depen- 
dence. Ignoring this energy-dependence allows us to 
write the solution of Eq. 116(1 for E < mx as 



n{E) 



n\{(TaV) 

m 



exp 



mx — E 



(21) 



For the interesting range of nix and E the exponential 
factor is always close to 1. Physically this represents the 
fact that positrons produced at energy E = mx will be 
lost from this "energy bin" primarily by down-scattering, 
not by annihilation. 

In order to predict the positron flux we use riH = 
1 cm~^. For the local dark matter density we use the 
canonical value pdm — 300 MeV cm^-^ = mx{nx + nx)- 
For the annihilation cross section we first consider an s- 
wave model with (aav) — (Tq where (Tq is fixed by the 
early- universe freeze-out calculation (Appendix . We 
compare the flux prediction with the best 95% CL upper 
limits in the 20-90 MeV range that are given in Fig. 4 of 
Ref. 28:]. At 20 MeV, the flux limit is approximately 
1.2 X IQ-^ cm-2 s-i sr-i MeV-\ for 30 MeV it is 
8.5 X 10"^ in these units, for 50 MeV it is 2.5 x 10"^ 
and for 90 MeV it is 9 x 10"'^. 

We find flux predictions exceeding these limits by fac- 
tors ^ 500. This implies that the annihilation cross 
section required in the early universe for reducing the 
particle density to the dark-matter level produces an un- 
acceptably large local positron flux and is thus excluded 
for mx > 20 MeV. For smaller mx we have no constraint. 
For a p-wave channel we do not obtain any limits because 
{aav"^) is strongly suppressed in the galaxy because of the 
small dark-matter velocities v of order 10"'^. 

Therefore, we confirm the conclusion of Refs. Q, |^ H, 
0,13 that an s-wave annihilation channel XX — > e+e~ is 
not acceptable for thermal relics. However, our argument 
does not depend on the uncertain dark-matter profile of 
the galactic bulge. 



IV. ENERGY TRANSFER IN STARS 

New MeV-mass particles could play an important role 
in stars. For example, they would be thermally produced 
in the collapsed core of a supernova and could contribute 
to the energy loss or the transfer of energy in these sys- 
tems [23. However, the AT-particles discussed here have 
stronger-than-weak interactions, implying that in a su- 
pernova these effects would be small compared to those 
of ordinary neutrinos. Therefore, even though the new 



particles would be thermally excited in a supernova core, 
there are no obvious observational consequences. 

In ordinary stars, and especially in our Sun, dark- 
matter particles will be trapped and contribute to the 
transfer of energy in potentially observable ways (3fl] . In 
the following we investigate if MeV-mass particles could 
be relevant in this context. The result is that for MeV- 
range masses the evaporation time is very short so that 
the steady-state abundance of X-particles in the Sun is 
too small to be important. 

A simple estimate of the energy conduction by the new 
particles can be worked out in a one-zone model of the 
Sun |3C|. One assumes that the dark-matter stationary 
distribution in the Sun is globally Maxwellian at a uni- 
form temperature Tx, with its maximum density found 
at a scale radius rx ■ Assuming that Tx is identical with 
the temperature at the solar center, and taking rx to be 
of order the solar radius, the luminosity carried by the 
new particles is of order 



10-*^ 



X 



^0 {as/ph)y/mx/MeV ' 



(22) 



where Nx is the total number of A-particles trapped in 
the Sun and as is the scattering cross section on electrons, 
taken to be comparable to the annihilation cross section 
as described in Appendix ^ For particles even more 
weakly interacting, i.e. for as ^ few pb, one enters the 
Knudsen regime where the effect of energy transfer is 
much smaller |^ [s^l ■ 

The steady-state number of dark-matter particles col- 
lected by the Sun arises from an equilibrium between 
capture and evaporation, i.e. Nx = A/Pe with A the 
number of particles captured per unit time and Pe the 
escape probability per unit time. We estimate the cap- 
ture rate to be ISS, ISJl 



A ~ 4.34 



nxG^MQRr, 



(23) 



where fgai ~ 300 km s""'^ is the mean square velocity of 
the galactic dark matter near the orbit of the Sun and 
nx its number density. 

The escape probability Pg is physically the ratio of the 
fraction of particles in the "escape-tail" of their distribu- 
tion to the typical time needed to repopulate it. A simple 
estimate is ^3 



GnMq 



2Tx 



mxvj 
' 2Tx 



(24) 



where Vf is the escape velocity from a typical posi- 
tion in the Sun. Near the solar center one has Vf 
(15.8Tx/TOp)^/^ m. Further, Sv ~ Vsms/mx where 
ms and Vs refer to the scattered particles, electrons in our 
case, for which we assume a thermal velocity distribution 
so that Vs = {'iTs/TnsY^'^ and we assume Tg = Tx- 

Based on these simple estimates and using a simplified 
solar model we find Nx ~ 10^^ (pb/cs)- Comparing this 
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result with Eq. (|22|l we conclude that the effect of the 
new particles is too small to be significant for the Sun. Of 
course, our estimate is rather crude considering that the 
escape lifetime of the dark matter particles is comparable 
to their orbital period. 



V. SUMMARY AND CONCLUSIONS 

We have analyzed some astrophysical and cosmological 
consequences of the intriguing possibility that the cosmic 
dark matter consists of MeV-mass particles [llllllll^. 
These particles are assumed to be thermal relics and thus 
have interaction cross sections that are larger than weak. 

Such particles do not have any apparent consequences 
for stellar evolution. In supernovae, their interaction is 
too strong so that neutrinos continue to play the lead- 
ing role for energy loss and energy transfer. In the Sun, 
the particle mass is so small that evaporation prevents 
the trapped steady-state population from growing large 
enough for a significant contribution to energy conduc- 
tion. 

We have derived new constraints coming from the 
low-energy positron component of cosmic rays. An s- 
wave annihilation cross section for XX —t e~^e~ as large 
as implied by the early-universe freeze-out calculation 
causes an excessive positron flux for mx = 20-90 MeV 
where experimental upper limits are available. There- 
fore, only p-wave annihilation is compatible with these 
constraints. These conclusions agree with those reached 

m Refs. laaaaii , but in our case they do not depend 

on the assumed dark-matter profile of the galactic bulge. 

The main subject of our work, however, was the im- 
pact of MeV-mass particles on big-bang nucleosynthesis 
(BBN). Significant modifications arise only for mx ^ 
20 MeV. The effects found are largely independent of the 
energy dependence of the annihilation cross section and 
even of its exact value, provided that it is high enough to 
take the new particles to equilibrium with the neutrinos 
or the electromagnetic plasma. 

For the neutrino-coupled case, the impact of the new 
particles is comparable to that of additional neutrino 
species. Notably, the primordial helium abundance is 
increased, exacerbating the discrepancy between predic- 
tions and observations. Therefore, such JT-particles are 
disfavored by BBN for masses up to about 10 MeV. 

For the electromagnetically coupled case, the BBN 
concordance would be severely disturbed for mx ^ 
2 MeV. However, there is a region mx = 4-10 MeV 
where the primordial helium mass fraction Yp is actu- 
ally reduced relative to the standard case while the pre- 
dicted remains compatible with observations. This 
non-trivial phenomenon is a consequence of the "entropy 
effect" discussed in the paper. While this effect was al- 
ready found in Ref. 'l3'|, it now assumes greater impor- 
tance because it slightly improves the discrepancy be- 
tween the BBN predictions for Yp and its observed value. 

In summary, the MeV-mass dark matter particles pro- 



posed in Refs. liliaSS 

as a source for positrons in the 
galactic bulge are not incompatible with BBN, provided 
the particle mass exceeds a few MeV. On the contrary, in 
the mass range mx — 4-10 MeV these particles slightly 
improve the concordance between BBN calculations and 
the observed helium abundance. It is quite fascinating 
that such exotic particles, far from being excluded, seem 
to have several beneficial consequences in astrophysics 
and cosmology. 
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APPENDIX A: RELIC DENSITY AND 
ANNIHILATION CROSS SECTION 



Once the mass mx and the spin multiplicity gx of 
a thermal relic are fixed, the requirement that it is the 
dark matter uniquely determines the annihilation freeze- 
out temperature Tp and the annihilation reaction rate at 
this temperature {(Jav)T-e- The relic mass density is 

^2 _ Sirh'^mxnx _ Sirrnxs (nx\ ...^ 
~ 3ml,m ~ 3{mpiH/hy \~ ) ' ^^^> 

where H/h = 100 km s^^ Mpc^^, mpi = 1.22 x 10^2 MeV 
is the Planck mass, and s is the entropy density. Note 
that for particles that are not self-conjugate fix only 
includes the particles so that f^DM = + ^x ■ 

Neglecting entropy-producing phenomena, the entropy 
per comoving volume remains constant so that 



\ s J Vs/Tf 



(A2) 



for T <Tp. If one assumes a functional form for {aav)T, 
an order-of-magnitude estimate of Tp is given by the con- 
dition nx{crv)TF = H{Tf). A more accurate formula 
is obtained by a semi-analytical treatment of the Boltz- 
mann equation [T^ 



where 



xo = log 



xp = xq — {n + i) logxo 



0.038 {n + 1) mpimxCTn 
J9*P 



Here, the parameterization 

{<Jav)T = cr„a;"" 



(A3) 



(A4) 



(A5) 
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with X = mx/T was assumed. The most interesting 
cases are obtained for n = (s-wave annihilation) and 
n = 1 (p-wave). The previous equations yield 



V s /o 
nxh^ 



3.79Vff;^(n + l)g;g+^ 
gstFnipimxCTn 



= 0.032 



gs*F(a-„/pb) 



(A6) 



with an accuracy better than 5%. Since H.uMh^ ~ 0.11 
[TtI. and xf — 15-20, it is easily seen that, for X ^ X, 
one has ctq '--^ 5 pb and ai ~ 10^ pb. 



APPENDIX B: KINETIC FREEZE-OUT 

We estimate the kinetic freeze-out temperature for X- 
particles coupled to the electromagnetic plasma. We 
make the rough approximation that the annihilation and 
scattering cross sections are comparable, da ^ (Jg, so that 
the ratio between the annihilation and scattering rates is 
essentially given by Ta/T^ ~ nx/ris- For our case the 
scattering targets are electrons and positrons so that 



\2ttz 



3/2 



^ cosh 



(Bl) 



where = fJ-e/T is the electron degeneracy parameter 
and z = nie/T. Therefore, as long as the X-particles are 
in kinetic equilibrium we have 



Fa _ 1 / mx_ 

Ts 2 I TOe 



3/2 



exp 



mx 



T 



1 



2 cosh 



(B2) 



This result applies at Tp only if {Ta/Ts)TF <C 1. As- 
suming mx = 1 MeV, one finds from Eq. IjASjl that 
Tp « 0.07 MeV and from a standard BEN code that 
^eiTp) « 0.32 X 10-^ so that (ra/r,)Tp « 4 X lO-'' which 
is indeed <C 1. 

Once a value for the annihilation effective cross sec- 
tion an has been determined, one can easily deduce a ki- 
netic freeze-out temperature Tk for the decoupling of the 
dark matter particles from the electromagnetic plasma, 
assuming that CTq ^ Gg. The condition Ts{Tyi) = H{Tk) 
implies 



g*{TK) 



mpins{T}<i) V 10.75 
With coshfe w 1 and « 3.36 one obtains 
ZK - (0.5 - e) logZK « 14.1 + log(CT,,^/pb) , 



(B3) 



(B4) 



where (agv) = Os^gz ^ . For typical values of Ug^n this 
gives Tk « 35 keV. 



APPENDIX C: ENTROPY GENERATION 
DURING DECOUPLING 

We sketch an argument that indicates that the entropy 
production associated with the dark-matter freeze-out is 
a sub- leading effect. To this end we assume that every 

e~ annihilation product is instantaneously thermalized 
and converted into photons. In this case the evolution 
equation for the X number density is |14| 



iix + SHnx 



-C[nx] 



where 



C[nx] = {(Tav) [n| - 71^ , 



(CI) 



(C2) 



The same equation applies to nx ~ nx because we al- 
ways assume equal distributions for both X and X. The 
subsequent e'^e^ annihilations imply 



+ 'SHn^ — +2C[nx] ■ 
The equilibrium energy densities are 



(C3) 



P-r 
Px 



2C{i)) 



4/3 



,4/3 



nx [mx + -T 



(C4) 



where the non-relativistic regime was used for the dark- 
matter particles. The energy injection by late annihila- 
tions is estimated as 



Pj 



T 



3 V30C(3) 
Px = ~mx 1 + 



C[nx] , 
3T 



2mx 



C[nx] 



(C5) 



The direct change of the total energy density due to this 
conversion of dark-matter particles into photons is a very 
small effect, as one can see by evaluating the ratios p^/p^ 
and hj/n^. 

The main non-negligible consequence is on the T-t rela- 
tion (first Friedmann equation) through the second term 
in the following equation 



dpc 



dt 



-3H{p + P), 



mx 



V45C(3) 



- 3 T 



C[nx] . (C6) 



The problem is then restricted to the solution of the equa- 
tion for nx in order to calculate the relevant quantity 
C[nx]- This was done by standard techniques after some 
standard substitutions (see e.g. Ref. 0|). We found a 
negligible effect on Fp, and a change of the and ^Li 
yields of order 0.1 in the units of Fig. [3 Therefore, we 
are indeed dealing with a sub- leading effect. 



10 



APPENDIX D: DEUTERIUM 
PHOTO-DISSOCIATION 



where U is the time at the onset of dissociation, tf is 
some late time where the D abundance is observed, and 



We show that the dissociation of fragile nuclides, prin- 
cipally deuterium, by late dark-matter annihilation is 
completely negligible. To this end we first follow stan- 
dard works on this subject (e.g. Ref. 35]) and note that 
the "first-generation" of nonthermal photons is rapidly 
degraded to have a spectrum with a high energy cut Eq ~ 
mg/22T because of the highly efficient pair creation and 
subsequent reactions on the background medium. This 
means that in order for some photons to survive and 
to be able to dissociate nuclei the temperature has to 
drop at least to values T < ml/22BD « 5 keV, where 
Bu Ri 2.2 MeV is the binding energy of deuterium. This 
implies that neglecting this phenomenon during the BBN 
epoch is certainly a good approximation. Note also that 
in our prefered models, MeV- mass X particles would not 
annihilate directly into photons, so the "first-generation" 
photon spectrum is already degraded in energy because it 
is produced by secondary effects of the primary electrons 
and positrons. 

Therefore, by simply looking at the expected energy 
spectrum, the most dangerous process would be the 
electro-disintegration of deuterium, for which a further 
0.5 MeV penalty in the energy should be considered as 
the electron rest mass would not be released. Without 
invoking a detailed analysis we conclude that our results 
up to mx ^2.7 MeV remain unchanged. 

For higher values of mx considered in our BBN anal- 
ysis, say 3-20 MeV, we reach the same conclusion by 
following the treatment of the late-time annihilations of 
a relic particle described in Ref. [s^. Assuming that 
electro-dissociation of deuterium is the only relevant phe- 
nomenon, the D depletion factor can be written as 



exp 



rEciT) 

TeD = S / dE fe {E, T) (TeD [E) . (D2) 



Here, aeD{E) « 13.3 ^lb[{E-Bl))/MeY]^/^ is the electro- 
dissociation cross-section of deuterium 's^. Further, 
fe {E, T) is the "steady state" dark-matter annihilation 
product spectrum. By arguments similar to the ones pre- 
sented in Ref. ^ol it is written as 



fe{E,T) « ar 



(. 



T 



(nxy 



\mx ) n-^Oe-^iE) \ s ' 
mxe[En,^AT) - E] 



(D3) 



(Dl) 



where aej{E) is the Compton scattering cross section 
and E^^y,{T) = Min[£;c(r), mx]- Note that the primary 
particle's "first" shower spectrum has roughly the same 
shape for 7 or e initiated showers [s^ [s^ . 

We have numerically solved the integral in Eq. ljDl|) 
by changing to the temperature variable and assuming 
a radiation dominated universe. The effects we found 
are completely negligible both for s- and p-wave annihi- 
lations. 

Note, however, that for showers induced by secondary 
photons one should replace in Eq. (jD3|l the huge fac- 
tor with rig, and the electro-dissociation cross section 
with the photo-dissociation one, while a penalty factor in 
the energy spectrum would enter. Making some extreme 
assumption one could thus gain up to a factor 10^^ of 
the previous estimates. Even so, it is not enough to af- 
fect by more than a few percent our simplified nuclides' 
predictions. The effects are even less pronounced for the 
preferred case n — 1 (p-wave annihilation). 
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